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We present calculations of the frequency-dependent spin susceptibility tensor of a two-dimensional electron 
gas with competing Rashba and Dresselhaus spin-orbit interaction. It is shown that the interplay between 
both types of spin-orbit coupling gives rise to an anisotropic spectral behavior of the spin density response 
function which is significantly different from that of vanishing Rashba or Dresselhaus case. Strong reso- 
nances are developed in the spin susceptibility as a consequence of the angular anisotropy of the energy 
spin-splitting. This characteristic optical modulable response may be useful to experimentally probe spin 
accumulation and spin density currents in such systems. 

Copyright line will be provided by the publisher 

Electrical manipulation of the electron and hole spins without the need of ferromagnetic materials and/or 
external magnetic fields is nowadays one of the central aspects in the field of spintronics. H] |2] [3] The 
presence of a sizeable spin-orbit interaction (SOI) in low-dimensional semiconductor structures and its 
modulation possibility (through electrical gating) make it a very prominent mechanism for the access and 
manipulation of the carriers spin states. 

It has been established that the dominant contributions to the SOI in quasi- two dimensional electron 
gases (2DEG) are the so called Rashba and Dresselhaus SO couplings. |4| The former results from the 
asymmetry of the confining potential that creates the 2DEG, while the latter arises due to the inversion 
asymmetry of the bulk. Several interesting effects and spin-based devices relying in these SOI mechanisms 
have been predicted and proposed in the last few years. For instance, the celebrated spin transistor proposed 
by Datta and Das |5|, and its recent non-ballistic version |6|. An intrinsic spin Hall effect in which a 
transverse spin current is driven by a dc electric field (without a net charge current) has been also predicted 
to occur due to SOI effects. |l7l[8][9 | More recently, a spin (Hall) accumulation has been observed through 
optical measurements ifTOl fTTl fT2]| . and lately, a purely electrical detection of a spin Hall current has been 
reported. |13| Electric-field-induced spin orientation in SOI coupled systems [141 [Tsl [161 [TtI [Tsll and 
strained semiconductors has been also explored. [ 19| 

On the other hand, the spin-splitting caused by SOI in electron systems opens the possibility of resonant 
effects via transitions between the spin-split states as a response to alternating electric fields. Il20ll2ni22l 
l23l l24l [25ll The importance of the study of such SOI effects in the dynamical regime (frequency depen- 
dent response) has been emphasized by several authors studying a variety of physical aspects, including 
spin and charge optical conductivities |22|, optical absorption spectra ||25ll26J , optical control of the spin 
Hall current through intense ac probing fields ||27l , electron-electron interaction effects |[28l[29l , electron- 
phonon interaction on spin Hall currents f30|, plasmon modes ||23l[3Tl[32l[33l , and the relation between 
the spin Hall conductivity and the spin susceptibility 1,28. .29. .341 or the dielectric function j|24ll . 
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The spin susceptibility plays a central role of the spin dynamics in a 2DEG. It gives the average spin 
polarization induced via electric-dipole or magnetic-dipole interactions. Thus, it can be used to obtain a 
magnetic susceptibility L28J or the electric-field-induced spin orientation factor. [ 16„ .35 1 Moreover, other 
transport properties like charge or spin Hall conductivities can also be expressed in terms of such spin 
response function. [34] 

Following S. I. Erlingsson et al.|34|, in this paper we report on the analytical and numerical calculations 
of the frequency-dependent spin susceptibility tensor of a 2DEG with Rashba and Dresselhaus SOI. In 
Ref. If34l expressions for the tensor components were obtained, however only approximated results were 
reported in the finite frequency regime. Their analytical expressions for the spin susceptibilities are valid 
as long as kso/kp << 1 and a << (3, where kso and are the characteristic spin-orbit coupling and 
Fermi wave numbers, while a and /3 are the Rashba and Dresselhaus SOI strength parameters, respectively. 

Here we show that in the more general case, particularly when there is a strong interplay between 
the Rashba and Dresselhaus SOI, very distinctive features of the optical spin susceptibility spectra arises 
in the system. This suggests that an optically modulable spin density response may be achievable in such 
systems. Furthermore, the calculated spectra show that the combination of the excitation at finite frequency 
and the interplay between the Rashba and Dresselhaus couplings could also be used for measuring the ratio 
between the SO coupling parameters. 

We consider a 2D free electron system lying at the z = plane subjected to spin-orbit interaction, with 
a Hamiltonian given by 

H= — ^-^-^^^-^1 + aik^ay - kya^) + (3{k^ax - kyffy) , (1) 

where k^ , ky are the components of the 2D electron wave vector, I is the 2x2 unit matrix and cr^ are 
the Pauli matrices. The second term corresponds to the Rashba SO coupling which originates from any 
source of structural inversion asymmetry (SIA) of the confining potential. The third term is the linear Dres- 
selhaus coupling which results from bulk-induced inversion asymmetry (BIA) in a naiTow [001] quantum 
well. Spin-orbit interaction appears as a relativistic coiTection (derived from the Dirac equation) to the 
Hamiltonian of a slow electron. It includes the gradient of a potential in which the electron moves. In 
atomic physics such term leads to the well known L • S coupling between the orbital and intrinsic angu- 
lar momentum due to the Coulomb potential. For an electron in a crystal environment there are several 
sources of potential gradient (impurities, confinement, boundaries, external electric field) which lead to an 
enhancement of SO coupling in solids. For quasi-2D systems the more significant contributions are those 
due to SIA (Rashba) and BIA (Dresselhaus). gl 

The eigenstates |kA) for the in-plane motion are specified by the wave vector k = {kx, ky) = k{cos9,s'm9) 
and chirality A = ± of the spin branches. The double sign corresponds to the upper {+) and lower (— ) 
parts of the energy spectrum given by 

where kso{d) = m*A{9)/h'^ is the characteristic SO momentum, A'^{9) = a'^ + (3^ — 2aPsm29 de- 
scribes the angular anisotropy of the spin splitting. At zero temperature, the two spin-split subbands are 
filled up to the same (positive) Fermi energy level sf but with different Fermi wave vectors q\{9) = 

\J 2m* ep/^^ + k1^{9) — \kso{9), determined from the equations e\{qx{9),9) = ep- Here, ep = 
?i^(fcQ — 2q^^)/2m* with /cq = \/27m being the Fermi wave vector of a spin-degenerate 2DEG with 
electron density n, and qso — m* \J + ffi jYi' . The SOI splits the Fermi line into two curves with radii 
given by q\{9) which, as the energy surfaces £A(k), are symmetric with respect to the (1,1) and (-1,1) 
directions in k— space (Fig. 1). When a or /? is null, the dispersions are isotropic and the Fermi contours 
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Fig. 1 Fermi contours qx(6) and the constant-energy-difference curve Cr{i^) defined by e+(k) — £-(k) = fuj, 
shown for two values of the photon energy, hcJi > hu>2- Cr{io) is a rotated ellipse with semi-axis of lengths (major) 
ka{uj) = hLj/2\a — P\ and (minor) kb{u}) = hu}/2\a+ 13\ oriented along the (1, 1) and ( — 1, 1) directions respectively. 
The sample parameters used here are n = 5 x lO^^cm"^, a = 1.6 x 10~^ eVcm, /3 = 0.5a and m* = 0.055m. 



are concentric circles. If a = ±(3 the spin-splitting along the (±1,1) direction vanishes and the spin- 
split dispersion branches are two circles with the same radius and displaced from the origin (along (t1 J) 
direction). 

Within the linear response Kubo formalism the spin susceptibility is given by 

X^^'H = ^ rfie^^"+"')*([a^(0,^^'(0)]> , M,m' = ^,2/ (3) 

where the symbol (• • • ) Sa / d^k f{£\(k)){- ■ ■ ) indicates quantum and thermal averaging, /(e) is the 
Fermi distribution function, and ?/ ^ 0+. This is a spin-spin response function for a spatially homogeneous 
(in-plane) perturbation oscillating at frequency oj. 

In the limit of vanishing temperature, eq. (O takes the form 

Xmm'M = A fd'k{^\a,\+){+\a,,\-) '^^!2^72 ^'L^,^ ^ ■ ,12 ' (4) 

71-^7 [e+[k) - e-[k)\^ - [n{uj + iri)\'' 

the prime on the integral indicates that integration is restricted to the region between the Fermi contours, 

q+{9) <k< q^ie), for which e_(k) < ep < e+(k), (Fig. 1). 
Using the result 

1 

= -i+lf^f^]-) = --rj^[Sf,xiacos9 - Psm9) + df,y{asme - Pcos9)] (5) 
the susceptibility tensor becomes 

2^ n rl-iO) Jc"^ 



1 r.nd^^AO) 



where 



A{9) Ak^A^i9)-[hiLu + z7^)Y 



5mm' (^) = S [5 (a COS 9 - (3 sin 9)^ + Sf_,y {a sin 9 - 13 cos 9)"^] 

+ ( 1 — Sfj_f^i){a cos 9 — (3 sin 9) {a sin 9 ~ (3 cos ( 
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It can be shown that Xxx{^j->) = Xyyi^) Xxyi^^) ~ Xyxi^)- Note also that for f3 = 0, Xxyi^^) 

Xyx{(^) =0. 

We can write the susceptibiHty in the form x^^' = x'^^, + «x^^/. The real part is 



where Tiil^ — \ep — Szfz{q±{0), 0)\ — 2q±{d)A{d) and the static value is 



Xmm'(O) = y 



(8) 



vo — m* /nh^ is the density of states of a spin-degenerate 2DEG, and 8(x) is the unit step function, 
e{x) = lifa; > Oande(x) = Oifx < 0. 
For the imaginary part we have 

X';.,'iu^) = 7rl'^{~\a,\+){+W,,\-)5{e+{k)~s^{k)-hu;) (9) 

= ^Jd9^^l^e[hi,-hn+{e)]e[hn^{e)^ni,]. (lo) 

These equations express the fact that the only transitions allowed between spin-split subbands e\ due to 
photon absorption at energy huj are those for which hn+{9) < Tau < hil.^{9). That is, for a given uj only 
those angles satisfying this condition must be considered in the integral (fTOb . see Fig. 2c. This is different 
to the pure Rashba (or Dresselhaus) case, where the whole interval [0, 27r] contributes to the integral for 
each allowed photon energy. The non-isotropic spin-splitting originated by the simultaneous presence of 
both coupling strengths, forces the optical excitation to be k— selective. 

In Fig. 2 we show Xxxi^^) as obtained from eqs. (fTTi-lfTOli. the xy component behaves similarly. The 
result is remarkably different from that of the pure Rashba or Dresselhaus case, where the spin-splitting 
is isotropic in the momentum space. For example, if /3 = 0, a 7^ 0, then x'^^' (^) = XRSfj.fj.' only for 
2ak^ < fiLo < 2afc_, otherwise it vanishes, where xr = ?iw/16a^, with fc± = '?±(/3 — 0) being 
independent of angle 9; (see Fig. 3). Thus, in this case the width A5 of the spectrum is A£/j = Aeji 
(or Aeb \f a — 0, j3 ^ 0); £r = m'a^ jf? and eu — m* jh^ are the SO characteristic energy 
scales for the Rashba and Dresselhaus coupling. As was discussed in Ref. Il22l . A£jj can be about an 
order of magnitude smaller than the width of the spectrum shown in Fig. 2b. Assuming that a > /5 and 
(fcso(6')/fco)^ < 1, we have A£ = 4/3fco + Af^ + A£d (if a < /3 the first term changes to Aako). Thus, 
the absorption bandwidth could be manipulated by tuning the coupling strength a and/or through variations 
of the electron density n = lt^/2TT. This fact could also be used to determine the sign of a — /3. ll22l 

To understand the structure of the spectra of Fig. 2, we further note that, according to eq. (O, the 
minimum (maximum) photon energy hLO+ (huj-) required to induce optical transitions between the initial 
A = — and the final A = + subband corresponds to the excitation of an electron with wave vector lying 
on the g+ (g_) Fermi line at 9+ = 7r/4 or 5tt/A (9^ — 37r/4 or 77r/4), giving fiuj± = Ml±{9±) = 
2fco|aT/3| T2m*(a=F/3)^/?i^. The absorption edges in the spectrum of Fig. 2b correspond exactly to ?iti>±. 
The function x^^/ (ijj) can also be written as a line integral along the arcs of the resonant curve Cr(w) lying 
within the region enclosed by the Fermi lines q\{9); see Fig. 1. The peaks observed in Fig. 2b correspond 
to electronic excitations involving states with allowed wave vectors on Cr{oj) such that |Vk(£+ — e-)\ 
takes its minimum value. The first (second) peak is at a photon energy Tiuja (fiuJb) for which the major 
(minor) semi-axis of the ellipse Cr{oj) (Fig. 1) coincides with the Fermi line q^{9+) hence 
huJa = hn_{9+) = 2ko\a~f3\+2m*{a-P)^/h'^ andhuJt = nn+{9-) = 2ko\a+ (3\-2m* {a+ (3)^ /h^ . 
The spectrum of x'xxi^) looks very similar to the joint density of states for the spin-split bands £+. ll22l The 
unequal splitting at the Fermi level along the symmetry (1, 1) and (—1, 1) directions is thus responsible for 
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Fig. 2 (c) Angular region (shaded) in k— space available for direct transitions as a function of photon energy. Only 
the shaded region contribute to the optical absorption [eq. (TO)]. The energy boundaries are given by hQ,+ {9) = 
ef — £- {q+{6), 9) and Ml- (6) = e+(g- {6),6)—eF- (b) Imaginary and (a) Real part of the optical spin susceptibility 
Xxx(<jj), vo ~ m* /nh? . For the frequencies = S7+(7r/4), ~ f2-(7r/4), oij, = f2+(37r/4), uj- = f2_(37r/4), 
see the text. The sample parameters are the same as in Fig. 1 . 



the peaks at photon energies TujJa and fiuj}, respectively, giving meaning to the structure of the spectrum. 
The overall magnitude and the asymmetric shape of the spectrum are due to the factor g^^i{6) / A'^{0) in 
eg. ([Toll. The results for several values of j3/a are shown in Fig. 3. 

The real part of XtJ-tJ.' {'^) presents additional spectral features. For photon energies in the range hu)a < 
Uuj < Tiwb we find numerically that it takes the constant values Xxa:(^) = ^'o/2andXxi/('^) = — (^'o/2)(a^ + 
0^)/2a(3. The spectral characteristics of the response displayed in Fig. 2a shows that the magnitude and 
the direction of the dynamic spin magnetization could be modified via electrical gating and/or by adjusting 
the exciting frequency. This suggests new possibilities of electrical manipulation of the spin orientation in 
a 2DEG in the presence of competing Rashba and Dresselhaus SO couplings. 

Following Ref. If36ll we have also obtained the static value of XtJ.ti' {'^) for finite momentum relaxation 
rate 77 > (see eq. ((Sjl). This parameter accounts phenomenologically for dissipation effects due to 
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Fig. 3 Imaginary part of the spin susceptibility x^^^^^) for several values of the ratio /3/q. Other parameters are as 
those used in Fig. 1 . 



impurity scattering. We found that, to linear order in er^d I'hf], it vanishes as (a 7^ /? 7^ 0) 



It is also possible to relate the spin current response to the spin density response. The definition of the spin 
conductivity {uj) describing a z— polarized-spin current flowing in the i— direction as a response to the 
field E{uj)y involves the commutator [J'f (t) , jy (0)], where ji = evi and J'f — h{az,Vi}/A are the charge 
and spin current operators, respectively. Using the velocity operator v(k) = V-i^H/h = ?ik/rn*+x(/3(T2; + 
aay)/h — y{a(Tx + (3(Ty)/h, this commutator can be written in terms of the correlators [ai{t) , aj (0)], 
( i = X, y), which determines the spin susceptibility Thus, the following relations can be derived 



e/87r Va^-/3V ^^0/2 - f3^ J i^q/2 

(13) 



<^vy{^) _f 2a(3 \ Xxxii^) , fa^+f3^\ Xxyii^) 



yy 

v^l2 ' W-py iyo/2 

These expressions are formally equivalent to eqs. (39) and (40) of Ref. ||34l . This connection is very 
convenient because a spin polarization is more experimentally accesible than a spin current. 

In summary, we have calculated the finite frequency spin susceptibility tensor of a two-dimensional 
electron gas with competing Rashba and Dresselhaus spin-orbit interaction. We find that the angular 
anisotropy of the energy spin-splitting introduced by the interplay between both SO coupling strengths 
yields a finite-frequency response with spectral features that are significantly different from that of a pure 
Rashba (Dresselhaus) coupling case. As a consequence, an optically modulable spin density response is 
then achievable in such systems which may be useful for spintronics applications. 
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